Majorana Fermions in Equilibrium and Driven Cold Atom Quantum Wires 
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We introduce a new approach to create and detect Majorana fermions using optically trapped ID 
fermionic atoms. In our proposed setup, two internal states of the atoms couple via an optical Raman 
transition — simultaneously inducing an effective spin-orbit interaction and magnetic field — while a 
background molecular BEC cloud generates s-wave pairing for the atoms. The resulting cold atom 
quantum wire supports Majorana fermions at phase boundaries between topologically trivial and 
nontrivial regions, as well as 'Floquet Majorana fermions' when the system is periodically driven. 
We analyze experimental parameters, detection schemes, and various imperfections. 



Majorana fermions (MFs), which unlike ordinary 
fermions are their own antiparticles, are widely sought for 
their exotic exchange statistics and potential for topolog- 
ical quantum information processing. Various promising 
proposals exist for creating MFs as quasiparticles in 2D 
systems, such as quantum Hall states with filling factor 
5/2 pp, p-wave superconductors [2], topological insula- 
tor/superconductor interfaces [3 [4], and semiconductor 
heterostructures [5-8 . In addition, MFs can even emerge 
in ID quantum wires, such as the spinless p-wave super- 
conducting chain [9] which is effectively realized in semi- 
conductor wire/bulk superconductor hybrid structures 
with spin-orbit interaction and strong magnetic field [TOl - 
[T3] . Although there are many theoretical and experimen- 
tal efforts to search for MFs, their unambiguous detection 
remains an outstanding challenge. 

Significant advances in cold atom experiments have 
opened up a new era of studying many-body quantum 
systems. Cold atoms not only sidestep the issue of dis- 
order which often plagues solid-state systems, but also 
benefit from tunable microwave and optical control of 
the Hamiltonian. In particular, recent experiments have 
demonstrated synthetic magnetic fields by introducing 
a spatially dependent optical coupling between different 
internal states of the atom [TH US] ? which can be gen- 
eralized to create non-Abelian gauge fields with careful 
design of optical couplings [I6j[l7]. For example, Rashba 
spin-orbit interaction can be generated in an optically 
coupled tripod- level system [18], which can be used for 
generating MFs in 2D [H [20] . 

In this Letter, we propose to create and detect MFs 
using optically trapped ID fermionic atoms. We show 
that optical Raman transition with photon recoil can in- 
duce both an effective spin- orbit interaction and an ef- 
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FIG. 1: (Color online.) (a) Optically trapped fermionic atoms 
form a ID quantum wire inside a 3D molecular BEC. Two Ra- 
man beams propagate along ki and directions, respectively. 
The recoil momentum k\ —ki — 2kx is parallel to the quan- 
tum wire, (b) Raman coupling between two fermionic states 
and a± induces a 2k momentum change from photon recoil, 
(c) RF-induced atom-molecular conversion. 



fective magnetic field. Combined with s-wave pairing in- 
duced by the surrounding BEC of Feshbach molecules, 
the cold atom quantum wire supports MFs at the bound- 
aries between topologically trivial and non-trivial super- 
conducting regions [lOl [11] . In contrast to the earlier 2D 
cold-atom MF proposals that require sophisticated op- 
tical control like tilted optical lattices [21] or multiple 
laser beams p~8j [20] , our scheme simply uses the Raman 
transition with photon recoil to obtain spin-orbit interac- 
tion. Moreover, compared with the solid-state proposals 
[3j [lOl [11] , the cold atom quantum wire offers various ad- 
vantages such as tunability of parameters and, crucially, 
much better control over disorder. 

Theoretical Model. We consider a system of optically 
trapped ID fermionic atoms inside a 3D molecular BEC 
(Fig. [TJ. The Hamiltonian for the system reads 

H = a t ( £ p + V + Srf ^ a P w 
p 
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The fermionic atoms with momentum p have two relevant 
internal states, represented by spinor a p = (a p ^,a p ^) T . 
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The kinetic energy is e p = |^ and optical trapping po- 
tential is V. As shown in Figs. [TJa) and (b), two laser 
beams Raman couple the states a p _k,i and a p +k,t with 
coupling strength B = Ql ^ 2 ; where £) e is the optical de- 
tuning, ^1(2) are Rabi frequencies, and k\ — k^ — 2kx 
is the photon recoil momentum parallel to the quan- 
tum wire. The bulk BEC consists of Feshbach molecules 
(b ^ + ai) [22-24 with macroscopic occupation in 
the ground state (bo) = S. The interaction between 
the fermionic atoms and Feshbach molecules can be in- 
duced by an RF field with Rabi frequency g and detuning 
Srf • The effective pairing energy is A = for fermionic 
atoms [25] . 

We can recast the Hamiltonian into a more trans- 
parent form by applying a unitary operation that in- 
duces a spin-dependent Galilean transformation, U = 
e ikfx(al it a x , t -4 4 a X); )d^ where x is the coordinate along 
the quantum wire. Depending on the spin, the transfor- 
mation changes the momentum by ±fc, Udp+k^U^ = a p ^ 
and Ua p -k^U^ = a p ,i- The transformed Hamiltonian 
closely resembles the semiconducting wire model studied 
in [TQl [11] and reads 

H = ^2 a t ( £ p - V + u P a z + Ba x ) a p + ^Aa^a^^ + h.c. 

p 

(2) 

where \i = — (Srf + V + £&) is the local chemical poten- 
tial and the velocity u = k/m determines the strength of 
the effective spin-orbit interaction. 

Topological and Trivial Phases. The physics of the 
quantum wire is determined by four parameters: the s- 
wave pairing energy A, the effective magnetic field £?, 
the chemical potential /i, and the spin-orbit interaction 
energy E so = mu 2 /2. For p ^ 0, the determinant of 
H p is positive definite, so the quantum wire system has 
an energy gap at non-zero momenta. For p = 0, how- 
ever, Hp yields an energy Eq = B — A 2 + \i 2 which 
vanishes when the quantity C = A 2 + pi? — B 2 equals 
zero, signaling a phase transition [10] [11] (see Fig. [2Jd). 
When C > the quantum wire realizes a trivial super- 
conducting phase. For example, when B <C A, \i all en- 
ergy gaps are dominated by the pairing term, yielding 
an ordinary spinful ID superconductor. When C < a 
topological superconducting state emerges. For instance, 
when B ^> A, /i, E so the physics is dominated by a single 
spin component with an effective p-wave pairing energy 
A p « A^; this is essentially Kitaev's spinless p-wave 
superconducting chain, which is topologically non-trivial 
and supports MFs [9]. 

With spatially dependent parameters (/i, B or A), we 
can create boundaries between topological and trivial 
phases. MFs will emerge at these boundaries [TDJ [TT] . 
Spatial dependence of \i (x) can be generated by addi- 
tional laser beams with non-uniform optical trapping po- 
tential V (x). Then C (x) can take positive or negative 
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FIG. 2: (Color online.) (a) Energy dispersion for spin-orbit- 
coupled fermions in a magnetic field. There is an avoided 
crossing at p = with energy splitting 2B (dark solid line). 
The horizontal dotted line represents \J A 2 + /i 2 , which has 
two crossing points when ^A 2 + p 2 < B (blue dotted line) 
and four crossing points when y/ A 2 + p 2 > B (orange dotted 
line), (b) Phase diagram for topological and trivial phases 
with respect to parameters of A and p. (c,d) C (x) can take 
positive or negative values, which divides the quantum wire 
into alternating regions of topological and trivial phases. 

' values, which divides the quantum wire into alternat- 
ing regions of topological and trivial phases [Figs. [2|c) 
and (d)]. Exactly one MF mode localizes at each phase 
boundary. The position of the MFs can be changed 
by adiabatically moving a blue-detuned laser beam that 
changes \i (x). Similarly, we can also use focused Raman 
beams to induce spatially dependent B (x) to control the 
locations of topological and trivial phases. 

Floquet MFs. It has been recently proposed that pe- 
riodically driven systems can host non-trivial topological 
orders [26] [27], which may even have unique behaviors 
with no analogue in static systems [28]. Our setup in- 
deed allows one to turn a trivial phase topological by 
introducing time dependence, generating 'Floquet MFs'. 
For concreteness we consider the time-dependent chemi- 
cal potential 

f/ix for t e [nT, (n + 1/2) T) 
M j "1^2 for t E [(n + 1/2) T, (n + 1) T) ' W 

which can be implemented by varying the optical trap po- 
tential V or the RF frequency detuning Srf- In addition, 
we assume the presence of a ID optical lattice which mod- 
ifies the kinetic energy e p — >• — 2 J cos (ka) cos(pa) and 
the spin-orbit interaction upa z — ^ 2 J sin (ka) sin (pa) a z 
in Eq.Q, where J is the tunnel matrix element and a is 
the lattice spacing. 

Let Hj be the Hamiltonian with \i = fij. The time- 
evolution operator after one period is then given by 
U T = e - iH * T l 2 e- iH ^ T l 2 . We define an effective Hamil- 
tonian from the relation Ut = e~ lHeiiT , and study the 
emergence of MFs in H e R. Eigenstates of i^ e ff are called 
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FIG. 3: (Color online.) Floquet MFs with two distinct flavors. 
Quasi-energy spectrum of H e ff and topological charges (Qo 
and Q^) are plotted for varying period T of the drive. Since 
the quasi-energy is defined up to an integer multiple of 2tt/T, 
it can support Floquet MFs at E = tt/T (thick red line) as 
well as E — (thick blue line) . The appearance of the two MF 
flavors is not necessarily correlated, and a single Floquet MF 
is present in much of the parameter space. The parameters 
are /ii = —J, /i 2 = — 3J, B = J, A = 2J, and 2ka — tt/4. 



Floquet states and represent stationary states of one pe- 
riod of evolution. The eigenvalues of H e R are called quasi- 
energies because they are only defined up to an integer 
multiple of 2tt/T. This feature, combined with the built- 
in particle-hole symmetry enjoyed by the Bogoliubov-de 
Gennes Hamiltonian, allows for Floquet MFs carrying 
non-zero quasi-energy. That is, since states with quasi- 
energy E and — E are related by particle-hole symmetry, 
states with E = or E — tt/T = —tt/T can be their own 
particle- hole conjugates. 

The existence of Floquet MFs is most easily revealed 
by plotting the quasi-energy spectrum of H e ^ in a finite 
system, which in practice can be created by introducing 
a confinement along the quantum wire. In Fig. |3j we 
plot the spectrum for a 100-site system with \i\ = — J, 
\±2 = —3 J, B = J, A = 2 J, 2ka = tt/4 for varying drive 
period T. Note that both H\ and H2 correspond to the 
trivial phase with 61,62 > 0. For small T, states with 
quasi-energy E = or E = tt/T are clearly absent from 
the spectrum — i.e., there are no Floquet MFs here. 

As one increases T, the gap at tt/T closes, and for 
larger T a single Floquet state with E = tt/T remains. 
We have numerically checked that the amplitude for this 
Floquet state peaks near the ends of the ID system. Thus 
it arises from two localized Floquet MFs and this state 
is associated with non-trivial topological charge Q n as 
we will see below. As one further increases T, another 
state at quasi-energy E = appears whose wavefunction 
again peaks near the two ends - a second type of Floquet 



MF - associated with a different, non-trivial topological 
charge Qo- Interestingly, the two flavors of Floquet MFs 
at E — and E = tt/T are separated in quasi-energies, 
and therefore, they are stable Floquet MFs as long as 
the periodicity of the drive is preserved. The presence of 
two particle-hole symmetric gaps changes the topological 
classification of the system from Z2 to Z2 x Z2. 

Two topological charges Qo and Q n are defined as 
follows. For translationally invariant quantum wire, 
the evolution operator has momentum decomposition 
Ut (t) = Y\ p Ut, p (t) for intermediate time r G [0,T]. 
After one evolution period, we have Ut = Ut (T) and 
Ut, p = Ut jP (T). The topological charge Qo ( or Qtt) is 
the parity of the total number of times that the eigenval- 
ues of Ut,o(t) and Ut^(t) cross 1 (or —1). The topolog- 
ical charges have the closed form 



?oQ 7r = Pf[M ] Pf[M w ] 
Qo = Pf[iVo] PfK], 



(4) 



where M p = log [Ut, p ] and N p = log [y / t/r^] are skew 
symmetric matrices associated with the evolution, and 
Pf[X] is the Pfaffian of matrix X. Here y/Uj\k is de- 
termined by the analytic continuation from the history 
of UT,k{r). Note that the product of topological charges 
QoQtt is analogous to the Z2 invariant suggested for static 
MFs [9] . In Fig. [3j we plot the topological charges Qo 
and Q n for various driving period T. Indeed, Floquet 
states at E = and E = tt/T appear in the range of T 
at which Qo and Q n equal to —1, respectively. 

Probing MFs. RF spectroscopy can be used to probe 
MFs in cold atom quantum wires [29l 32 . In particu- 
lar, we consider spatially resolved RF spectroscopy [33] 
as an analog of the STM. The idea is to use another 
probe RF field to induce a single particle excitation from 
the fermionic state (say a a ) to an unoccupied fluores- 
cent probe state /. Contrary to conventional RF spec- 
troscopy, a tightly confined optical lattice strongly local- 
izes the atomic state /, yielding a flat energy band for 
this state. By imaging the population in state /, we gain 
new spatial information about the local density of states. 

For example, by applying a weak probe RF field 
with detuning S RF from the a a -f transition, the pop- 
ulation change in state / can be computed from the 
linear response theory I (x,v) = ^ (f^ (x) f (#)) oc 
pa a (x, —Ji (x) — 5' RF + s)Q (ft (x) + 5' RF — e). Since the 
MFs have zero energy in the band gap and are spatially 
localized at the end of the quantum wire, there will be an 
enhanced population transfer to state / with frequency 
5' RF = £ — \i (x*) at the phase boundary x*. If the a a - 
f transition has good coherence, we can use a resonant 
RF 7r-pulse to transfer the zero-energy population from 
a a to /, and then use ionization or in situ imaging tech- 
niques [3U |35] to reliably readout the population in / 
with single particle resolution. Floquet MFs can also be 
detected in a similar fashion. Since a Floquet state at 
quasi-energy E is the superposition of energy states with 
energies E + 2utt/T for integer n, we should find the 
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Floquet MFs at energies (or ir) + 2nir /T for (or n) 
quasi-energy Floquet MFs, respectively. 

Parameters and Imperfections. We now estimate the 
experimental parameters for cold atom quantum wires. 
(1) The spin-orbit interaction energy is E so — mu 2 /2 < 
^rec,o? with recoil energy E rec « 30 (2ir) kHz for 6 Li 
atoms. If we use n sequential A transitions, the spin-orbit 
interaction strength can be increased to = nk/m and 

Eg 7 ^ = n 2 E so . (2) The s-wave pairing energy A = gS 
can be comparable to the BEC transition temperature 
kT c ~ h 2 n^ 3 /m before the BEC is locally depleted. For 
molecule density no = 10 cm" 3 [23|l24j. we have A ~ 
10 (27r)kHz. (3) The effective magnetic field B = 5^ 

and the depth of the optical trap Vb ~ -j- can be MHz, 
by choosing large detuning 5 ~ 100 (27r)THz and Rabi 
frequencies Q ~ 50 (27r)GHz, while still maintaining a 
low optical scattering rate r « ^-7 ~ 1 (27r)Hz. (4) The 
transverse oscillation frequency of the ID optical trap can 

be v ~ \J~^2 ~ 150 (27r)kHz for a laser beam with waist 
w = 15/im. Since v is much larger than the energy scales 
of E so and |A|, it is a good approximation to consider a 
single transverse mode of the quantum wire. 

In practice, there are various imperfections such as par- 
ticle losses due to collision and photon scattering, finite 
temperature of BEC, and multiple transverse modes of 
the quantum wire. (1) The lifetime associated with pho- 
ton scattering induced loss can be improved to seconds 
using large detuning and strong laser intensity, and the 
collision-induced loss can be suppressed by adding a ID 
optical lattice to the quantum wire. (2) At finite tem- 
perature the BEC order parameter will fluctuate, and 



the effects can be examined by considering a spatially 
dependent order parameter Eoe 1 ^^. A large phase gra- 
dient (j) x = d<j)/dx > 2 ^ will close the energy gap, and 
the MFs will merge into the continuum. To sustain the 
energy gap, the fluctuations in the phase gradient should 
be small, i.e., yflp x (T)> 

thermal 

< 4W X (T*)> 

thermal 

^=^, with critical temperature T*. Thus, the BEC tem- 
perature should be below min {T*, T c } ~ 50nK. (3) Since 
the quantum wire has a finite transverse confinement, 
other transverse modes might be occupied and coupled 
non-resonant ly. Nevertheless, recent numerical and ana- 
lytical studies [El [131 |36l [37] show that MFs can be ro- 
bust even in the presence of multiple transverse modes, as 
long as an odd number of transverse quantization chan- 
nels are occupied. These results may potentially relax the 
requirement of tight confinement of the quantum wire. 

In conclusion, we have proposed a scheme to create and 
probe MFs in cold atom quantum wires, and suggested 
the creation of two non-degenerate flavors of Floquet MF 
at a single edge. We estimated the experimental param- 
eters to realize such implementation, considered schemes 
to probe for MFs, and analyzed imperfections from re- 
alistic considerations. Recently, it has been discovered 
that braiding of non-Abelian anyons can be achieved in 
networks of ID quantum wires [38 , which would be very 
interesting to explore in the cold atoms context. 
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